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Abstract 

I have analyzed reduced neutron widths (r°) for the subset of 1245 resonances in the nuclear 
data ensemble (NDE) for which they have been reported. Random matrix theory (RMT) predicts 
for the Gaussian orthogonal ensemble (GOE) that these widths should follow a x 2 distribution 
having one degree of freedom (v = 1) - the Porter Thomas distribution (PTD). Careful analysis of 
the values in the NDE rejects the validity of the PTD with a statistical significance of at least 
99.97% {y = 0.801 ±0.052). This striking disagreement with the RMT prediction is most likely due 
to the inclusion of significant p-wave contamination to the supposedly pure s-wave NDE. When an 
energy dependent threshold is used to remove the p-wave contamination, the PTD is still rejected 
with a statistical significance of at least 98.17% (y = 1.217 ± 0.092). Furthermore, examination of 
the primary references for the NDE reveals that many resonances in most of the individual data 
sets were selected using methods derived from RMT. Therefore, using the full NDE data set to 
test RMT predictions seems highly questionable. These results cast very serious doubt on claims 
that the NDE represents a striking confirmation of RMT. 



*Electronic address: 



koehlerpe@ornl.gov 



I. INTRODUCTION 

The nuclear data ensemble (NDE) is a set of 1407 resonance energies consisting 

of 30 sequences in 27 different nuclides. The ensemble was assembled to test predictions of 
random matrix theory (RMT) [3J. Fluctuation properties of resonance energies in the NDE 
were found to be in remarkably close agreement with RMT predictions for the Gaussian 
orthogonal ensemble (GOE). Although there have been several other successful tests of RMT 
using nuclear resonances (e.g., {4, 5]), the NDE is perhaps the most important because, as 
stated in Ref. {3]], "As a result of these analyses, it became generally accepted that proton 
and neutron resonances in medium weight and heavy nuclei agree with GOE predictions." 
Hence, the NDE routinely is cited as providing striking confirmation of RMT. 

Reduced neutron widths (r°) have been reported for a subset of 1245 resonances in the 
NDE of Ref. {2], consisting of 14 to 178 measurements for 24 nuclides. If the GOE correctly 
describes the data, RMT predicts these widths should follow a \ 2 distribution having one 



degree of freedom (z/ = 1) - the 



3 orter Thomas distribution (PTD). It has be argued that 

1 n 

6] and more robust |7J than other RMT predictions. Also, 



the PTD is more generally valid 

it has been demonstrated js] that a much more reliable analysis of spectra fluctuations can 
be performed using neutron widths than using resonance energies. In addition, it has been 
shown {9] that resonance width fluctuations are more sensitive than energy fluctuations to the 
degree of chaos in model quantum systems. In addition, it is straightforward [lOj to account 
for experimental effects such as missed or spurious resonances, and to use the statistically 
efficient maximum-likelihood (ML) method while using width data, but the same cannot be 
said for energy data. For these reasons, a test of the PTD using the NDE neutron widths 
could be very valuable. 

In Ref. j2| such a test on a subset of the NDE data was described from which it was 
concluded that there was satisfactory accord between theory and experiment. However, 
there are several problems with the analysis of Ref. j^, which I will describe below. I find 
that when the data are analyzed more carefully, they do not agree with the PTD. 

I will show below that the NDE suffers from significant p-wave contamination. For such 
resonances, reported reduced neutron widths T° are actually only "effective" reduced widths 
gY n j\[E^, where g = ^j+ih^+i) (with J, /, and j being spins of the resonance, target, and 



(2/+l)(2j+l) 

neutron, respectively) is the spin statistical factor. All s-wave resonances for NDE nuclides 



2 



included in my analysis have J = ~ and hence g = 1 and g»r° = r°. However, there are two 
resonance spins possible for p-wave resonances for these nuclides; J — | or |, and hence 
g = 1 or 2. As a reminder of these facts, I will use gT° n when referring to reduced neutrons 
widths in the remainder of this paper. 



II. RECONSTRUCTING THE NDE 



To my knowledge, the resonance energies and neutron widths in the NDE have never 
been published as a set. It should be possible, however, to reconstruct this information from 

nn 

the nuclides and corresponding number of resonances given in the NDE papers 121] and the 
primary references cited therein. Unfortunately, several of the references listed in the NDE 
papers are "private communication", there are a few errors in the number of resonances 
reported in Refs. and this reference does not explain why some of the nuclides in their 
NDE were excluded from their analysis of the neutron widths. For example, in Ref. Q it 
is stated that 1182 neutron widths in 21 nuclides were included in the analysis. However, 
the sum of the number of resonances in these 21 nuclides reported in this same reference is 
actually 1194. One resonance in 182 W was reported without a neutron width in the primary 
reference [ll|, so only 40 of the 41 resonances for this nuclide can be included in the width 
analysis. Therefore, the total number of widths is actually 1193, not 1182. Also, 19, 54, 47, 
and 21 resonances have been reported for 154,156,158,160 Gd; respectively) not 1 9; 47> 2 1, 

and 54, respectively as stated in Ref. {2]. In addition, there are three nuclides ( 160 Dy, 164 Dy, 



and 186 W, with 18, 20, and 14 resonances, respectively) in the NDE of Ref. l2| that were 

nn 

not included in their width analysis, even though neutron widths were available [TJJ, [13J . 

In total then, the NDE of Ref. [2] should contain 1245 neutron resonances in 24 nuclides 
for which widths have been reported. I obtained resonance energies and neutron widths for 
these nuclides from the primary references given in Table HI and cross checked these data 
with information in the EXFOR/CSISRS 14j database. Given the problems noted above, 
I cannot be certain that I have analyzed the same data as those in the NDE of Ref. j^J. 

lowever, using data from the primary references should minimize any differences with Ref. 

2| as well as make it easier for others to reconstruct the data set I have used. 
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TABLE I: NDE nuclides. 



Nuclide Ref. E max (keV) ML Results 

Minimum Threshold p-free Threshold 

















T 

- 1 max 






«A r-M 

64 Zn 


[15, JJ 


367.55 


0.024 


103 


1. 


.35 


+0.24 
—0.22 


0.05 


99 


1.54 


+0.29 
—0.26 


bb Zn 


[15, iH 


297.63 


0.025 


65 


0. 


.68 


+0.23 
—0.21 


0.05 


61 


0.74 


+0.27 
—0.25 


b8 Zn 


[15,1J 


247.20 


0.019 


45 


0. 


.75 


+0.27 
—0.25 


0.05 


41 


0.95 


+0.36 
—0.32 


114 Cd 


[12] 


3.3336 


0.137 


17 


0. 


.35 


+0.54 
—0.34 


0.45 


11 


2.0 


+1.5 
— 1.2 


152 Sm 


[18] 


3.665 


0.025 


70 


1. 


.14 


+0.27 
—0.25 





.1 


62 


1.55 


+0.40 
-0.38 


154 Sm 


[18] 


3.0468 


0.036 


27 


0. 


.76 


+0.38 
—0.36 





.1 


22 


1.32 


+0.65 
—0.55 


154 Gd 


[12] 


0.2692 


0.15 


19 


0. 


.44 


+0.58 

— U.4o 





.2 


18 


0.49 


+0.64 

n a q 
— U.4o 


156 Gd 


f8] 


1.9908 


0.009 


54 


1. 
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— U.zo 





.2 
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— U.4y 


1&8 Gd 
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47 


0. 


.75 


+0.25 
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.2 


35 


1.17 


+0.54 
—0.47 


ibU Gd 


[12] 


3.9316 


0.013 


21 


0. 


.55 


+0.34 
—0.33 





.2 


16 


0.83 


+0.75 
—0.65 




[13] 


0.4301 


0.07 


18 


0. 


.83 


+0.57 

— U.Ol 





.2 


14 


1.41 


+1.0 
— U.oo 


lb2 Dy 


[13] 


2.9572 


0.046 


46 


1. 


.02 


+0.33 
—0.32 





.2 


40 


0.99 


+0.47 
—0.43 


u>i Dy 


[13] 


2.9687 


0.04 


20 


0. 


.82 


+0.51 
—0.44 





.2 


16 


2.3 


+1.2 
— 1.0 


Er 


[19] 


4.1693 


0.02 


109 


0. 


.85 


+0.18 
—0.17 





.3 


78 


1.85 


+0.49 
—0.45 


Er 


[19] 


4.6711 


0.04 


48 


0. 


.80 


+0.30 
—0.27 





.3 


37 


1.32 


+0.62 
—0.55 


1 70-1-1 


[19] 


4.7151 


0.05 


31 


0. 


.36 


+0.34 
—0.31 


.03 


17 


3.6 


+1.6 
— 1.3 
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[20] 
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0.05 


55 
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0. 


06 
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-0.26 
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[20] 
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0.02 


19 


0. 


.80 
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0. 


06 


16 


1.29 
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i76 Yb 


[20] 
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0.015 


23 


0. 


.04 


+0.29 
-0.03 


0. 


06 


15 


1.05 
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-0.55 


182 W 


M 


2.6071 


0.069 


40 


0. 


.76 


+0.37 
-0.35 


0. 


15 


34 


1.50 


+0.62 
-0.55 


184 W 


M 


2.6208 


0.04 


30 


0. 


.62 


+0.34 
-0.31 


0. 


15 


26 


0.99 


+0.54 
-0.48 


186 W 




1.1871 


0.07 


14 


1. 


.23 


+0.78 
-0.62 


0. 


15 


13 


1.32 


+0.93 
-0.75 


232rpL 


[21] 


2.988 


0.016 
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0. 


.76 


+0.13 
-0.12 


0. 


26 


123 


1.78 


+0.36 
-0.34 


238TJ 


[21] 


3.0151 


0.0045 


146 


0.79 ± 0.12 


0. 


47 


84 


1.02 


+0.39 
-0.34 


W.A. 








1245 0.801 ± 0.052 






978 


1.217 ± 0.092 
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III. IMPORTANCE OF THRESHOLD IN ML ANALYSIS OF NEUTRON 
WIDTHS 



Because the PTD is a special case (degrees-of-freedom is = 1) of the family of x 2 dis- 
tributions, it is assumed that the data are distributed accordingly and the ML method is 
used to estimate the most likely value of is. Ideally, the data should be complete (no missing 
resonances) and pure (all resonances have the same parity). Unfortunately, all experiments 
from which neutron widths have been determined have finite thresholds for detecting reso- 
nances and for separating small s- from large p-w&ve resonances. Not properly accounting 
for these effects can result in substantial systematic errors in ML estimates of v. Given 
the shape of the x 2 distribution as a function of is, neglecting the effect of missed s-wave 
resonances below threshold will, in general, lead to a falsely large value of v from the ML 
analysis. Conversely, including even just a few p-wave resonances in an s-wave set will, in 
general, lead to a falsely small value of is. This potential problem is especially important 
for many of the NDE nuclides because they are near the peaks of the p- and minimum of 
the s-wave neutron strength functions, and therefore neutron width is a much less reliable 
indicator of resonance parity. 

Below I will show that the NDE is seriously contaminated by p-wave resonances. It is 
also almost certainly incomplete, as illustrated in Figs. [T]and[2j Reduced neutron widths for 
all 1245 resonances in the NDE are shown as a function of resonance energy in Fig. [H from 
which it can be seen that their are fewer small widths as the energy increases. This is just 
what is expected from well-known [22| instrumental effects that decrease sensitivity as the 
energy increases; hence, more small resonances are missed at higher energies. This is further 
illustrated in Fig. [2] in which distributions of reduced neutron widths are shown for the 100 
lowest- and highest-energy NDE resonances. As can be seen in this figure, the distribution 
of the highest-energy set is substantially narrower (and hence is in better agreement with 
a x 2 distribution having a larger is) than the lowest energy one. Again, this is just what is 
expected if more resonances are missed as energy increases. 

Typically, these difficulties have been surmounted by using a energy-independent thresh- 
old as an integral part of the ML analysis, implicitly assuming that all s-wave resonances 
above threshold have been observed. We recently have shown [23( that an energy- dependent 
threshold (on gT^) of the form T = aE n , where a is a constant factor, offers three ad- 
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vantages compared to using a energy-independent threshold. First, p-wave contamination 
is eliminated equally effectively at all energies. This is because the penetrability factor for 
p-waves differs from s-waves by (to good approximation) a factor of E n . Second, experiment 
thresholds have approximately this same energy dependence; thus, possible diffusiveness of 
the instrumental threshold can be surmounted equally effectively at all energies. Third, 
statistical precision of the analysis is maximized by allowing the largest j?-wave-free set of 
s-wave resonances to be included. 

Analyzing a data set comprised of many different nuclides such as the NDE involves 
at least two additional potential pitfalls. First, as is evident from Fig. [31 the apparent 
sensitivities of the various experiments from which the NDE was derived differ by several 
orders of magnitude. Therefore, if the entire NDE were analyzed as a single set (as was 
done in Ref. |2j), the threshold must be at least as high as the highest apparent individual 
threshold. But doing this will substantially reduce the statistical precision of the result 
and at least partially negate the reason for assembling the NDE in the first place. Second, 
for a given set of data and threshold, the ML-estimated average reduced neutron width 
might be substantially different from the one estimated from the data (e.g., from a simple 
average or by assuming v — 1). Therefore, it is important to include the average reduced 
width as a parameter in the ML analysis. Furthermore, the difference between these two 
estimated values will likely vary from one NDE nuclide to the next. For these two reasons, 
it is important that separate ML analyses be made for each NDE nuclide. 

To minimize the effects of the above problems, I have done separate ML analyses for 
each nuclide in the NDE using (a range of) separate energy-dependent thresholds, and then 
combined these individual results for comparison to theory. 



IV. AN IMPROVED ML ANALYSIS OF THE NDE NEUTRON WIDTHS 



The analysis technique was briefly described in Ref. 23|]. Each resonance A has an energy 
Ex and reduced neutron width gT^ n . The probability density function (PDF) f{x\i>) for a 
X 2 distribution is given by, 



f(x\u)dx 



2G(|) 



,y) 2 exp( 



— )ax, 
2 ' ' 



(1) 
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FIG. 1: (Color online) All 1245 reduced neutron widths in the NDE (red X's). Data for each 
nuclide have been normalized to their respective average reduced neutron widths and maximum 
energies. Blue circles depict those resonances which have been identified as being p wave or of 
uncertain parity. 




FIG. 2: (Color online) Reduced neutron width distributions for the 100 lowest- (blue cirles) and 
highest- (red X's) energy resonances in the NDE. The blue solid curve depicts the RMT prediction 
for the GOE (the PTD) whereas the red dashed curve is for a x 2 distribution having i/=1.5. See 
text for details. 
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FIG. 3: (Color online) Blue circles depict minimum reduced neutron widths, normalized to their 
respective averages, for each of the nuclides in the NDE, versus mass number. 

where G(|) is the gamma function for |, and x — > gr° n /E[gr°J, where Efgr^J is the 
expectation (average) value of the reduced neutron width, with E[«] denoting the expectation 
value operator. The ML method is used to estimate most likely values of v and E[F° n ] as 
well as their uncertainties. 

To facilitate comparison of the various nuclides in the NDE (e.g., FigsiU El El and[5j and 
Tabled]), I have normalized data for each nuclide to their resp ective average reduced widths 



< gT° n > (as reported in the primary references or in Ref. [24|]) an d maximum energies E max . 
Hence, thresholds (on gT^/ < gT° n >) can be expressed as T' = T max E n / E max and thus, 
T max is the maximum value of the threshold relative to the average reduced neutron width. 
Values of T max used in the analyses are given in Table fl] 

The joint PDF for statistical variables gT^ n and E\ is defined in a 2D region X given by 
inequalities E\ < -E max and gT° Xn > T(E\). The expression for this PDF reads 



^ (^^r An |z.,E[^ n ])=c/ 



V 



(2) 



The factor C, ensuring a unit norm of h°, is v- and E[gr° n ]-dependent. The ML function 
was calculated from all n pairs [E^ xp , gT^] , obtained from the experiment, which fall into 



8 



the region X. Specifically, 

no 

L MWJ) =Yl h °(K P i9 T °C P I ^E«J) . (3) 

i=l 

For the initial analyses, thresholds just below the smallest observed resonance for each 
nuclide were used. ML results with these thresholds are given in column z/ m ; n of Table [H 
Because experiment thresholds might not be precisely sharp, it is expected that the resulting 
z/ min values would be systematically a bit large. However, almost all z/ min values are less than 
the PTD value of 1.0 . 

Contour plots of ML functions for 232 Th, calculated at two different thresholds, in the 
form 

z(y, E[gT° Xn ]) = 2^ [l n L max -hxL (v, E[gT° Xn })] 1 (4) 

are shown in Fig.Hl Here, L max is the maximum of the ML function. Contours at fixed z = k 
encircle approximately the ka confidence region, and were used to derive the la uncertainties 
given in Table [H Careful statistical analysis in Ref. 23|] verified that these contours are 
reliable. The weighted average of results at minimum threshold for the 24 nuclides in the 
NDE is v — 0.801 ±0.052, which is 3.8 standard deviations smaller than the predicted result 
of v = 1. Hence, these data reject the PTD with a statistical significance of 99.98%. To 
check this result, the combined probability for the 24 NDE nuclides was calculated using 
Fisher's and Stouffer's (both weighted and unweighted) techniques [25[. These methods 
yielded combined confidence levels of 99.97%, 99.98%, and 99.99%, respectively, in good 
agreement with the weighted-average result. Hence, at minimum thresholds the NDE data 
reject the PTD with high confidence and indicate that v is significantly less than 1.0. h v 



value significantly less than the PTD could be a sign of interesting physics 



26 



-|29|. 



However, a more likely explanation is that the NDE contains sizeable p-wave contamination. 



V. CLEANSING THE NDE OF P- WAVES 



Results of the maximum-likelihood analyses for three NDE nuclides are shown in Fig. 
|5j On the left of this figure, normalized reduced neutron widths are plotted as functions of 
normalized resonance energy. On the right side of this figure, v values from the ML analyses 
are plotted versus threshold coefficients T max . 
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FIG. 4: Plots of z(u, E[gF° Xn \) constructed from ML analyses, at minimum (left) and p-wave-free 
(right) thresholds, of the NDE 232 Th data. On each plot, a filled circle indicates the location of 
L maj( , and a dashed horizontal line is drawn at v = 1, the PTD value. 



That the NDE is contaminated by £>-wave resonances is evident from this figure in two 
ways and from Figs. [1] and |6j First, resonances in the NDE that have been identified (in 



Rcfs. 



24 



30j and references contained therein) as p-wave or of uncertain parity are shown 



as open circles in Figs. [T]and[5j For example, in Ref. 31], 58 p-wave resonances in 232 Th 
were assigned on the basis of 7-cascade information, 13 of which are in the NDE. One of 
these 13 resonances also is known (32), 33] to be p wave by its observed parity- violating 
asymmetry. Percentages of p-wave or unknown-parity resonances for NDE nuclides, as a 
function of mass number, are shown in Fig. |6j For over half (14/24) of the NDE nuclides 
analyzed herein, these resonances account for 5% or more of the total, for 10 of the 24 they 
are at least 10%, and in the three worst cases about 35%. 
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FIG. 5: (Color online) Left: normalized reduced neutron widths versus normalized resonance 
energies for 114 Cd, 166 Er, and 232 Th resonances in the NDE. Red X's depict all resonances in the 
NDE whereas blue circles show resonances previously identified as being p wave or of uncertain 
parity. Right: Red X's depict v values from ML analyses versus thresholds used, for the same 
three nuclides. Error bars represent la confidence levels. Black dashed vertical lines correspond 
to thresholds depicted by black dashed curves in the left part of this figure. See text for details. 

Second, that many of the NDE resonances are in fact p-wave is reinforced by the behavior 
of the v values from the ML analyses as functions of threshold, as shown in the right side of 
Fig. In all three cases shown, v systematically increases with threshold before gradually 
stabilizing. This is just the behavior expected for a population of s- wave resonances contam- 
inated by p-wave resonances. Similar fractions of previously identified p-wave resonances 
and trends in v with threshold are seen for several of the other NDE nuclides. 
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FIG. 6: (Color online) Blue circles depict percentages of NDE resonances for each nuclide which 
are known to be p wave or are of uncertain parity, as a function of mass number. 

Removing effects of these p-w&ve resonances from the NDE ML analysis is a simple matter 
of raising thresholds until they are above the largest previously-identified p-w&ve resonance 
and/or v stabilizes as a function of threshold. Typical "p-wave free" thresholds for three 
NDE nuclides are shown as dashed curves in the left-hand part of Fig. and corresponding 
values of T max are depicted by dashed vertical lines in the right-hand part of this figure. 
Degrees-of- freedom values for each of the NDE nuclides at these "j?-wave free" thresholds 
are given in column five (v p f) of Tabled! The resulting weighted average for the NDE is still 
in conflict with the RMT prediction for the GOE, albeit in the opposite direction from the 
result using the lowest thresholds: v = 1.217 ± 0.092, corresponding to a confidence level of 
98.17% for excluding the PTD. Fisher's and unweigthed and weighted Stouffer's confidence 
levels are somewhat higher; 99.15%, 99.81%, and 99.37%, respectively. Hence, when the 
NDE is cleansed of p-wave resonances, the data still reject the PTD with high confidence. 



VI. DISCUSSION 



I have shown that when neutron widths in the NDE are analyzed carefully and in such a 
way as to eliminate p-wave resonances, the data exclude the PTD with fairly high confidence 
{y = 1.217 ± 0.092). At the same time, it has been shown [l|, [2] that the NDE as a whole 
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agrees remarkably well with energy-level fluctuations predicted by RMT for the GOE. Given 



the greater sensitivity of widths over energies to the degree o 



27 



29 



34 



chaos in the system [9|, and 
351 ] . perhaps this dichotomy can 



to effects related to the "openess" of the system 
be resolved. However, I am not aware of any published model which does so while at the 
same time yielding v > 1. 

For example, in the calculations of Ref. [9j, collective effects sometimes resulted in regions 
of model space where transition-strength distributions deviated strongly from the PTD but, 
at the same time, energy-level fluctuations were consistent with GOE predictions. How- 
ever, transition-strength distributions in this model deviated from the PTD in the direction 
opposite (v < 1) of the NDE. 

Reaction effects also should be considered to explain the dichotomy. For example, it 
recently was proposed [35( that the standard transformation of measured to reduced s- 
wave neutron widths, T° = T n / y/E, should be modified for nuclides near the peaks of 
the s-wave neutron strength function. This newly proposed transformation, T° = T n / y/E x 
7r(^) 1 / 2 (_E+|i?o|), where E is the energy (relative to threshold) of the s-wave single-particle 
state responsible for the peak in the neutron strength function, could affect (relative to the 
standard transformation) the shape of the T° distribution for values of E near the resolved- 
resonance energy range. However, as most of the NDE nuclides are not near the peaks of the 
s-wave neutron strength function, this proposal should not affect the NDE, and therefore 
cannot explain the dichotomy. 



Several explanations related to the "openess" of the system [271, |28|, |34| and to correla- 



tions between incoming and outgoing channels 29] have been proposed for the observation 



231 ] that T° distributions for 192 > 194 pt resonances are significantly broader [y 0.5) than 



the PTD. However, width distributions 



2a 



29j predicted by these theories are better char- 



acterized by v < 1 rather than v > 1 as observed for the NDE. 

Considering how data in the NDE were selected suggests another solution. It often has 



been stated that A3 is very sensitive to missing or misassigned resonances (see e.g. 



36- 



for recent work on this subject). Therefore, if the NDE was pure but incomplete, or 
vice versa, it would not be expected to agree well with the spacing statistics used in Refs. 

. However, it seems plausible that because the NDE is both impure and incomplete, 
it can be made to agree with these statistics, especially considering that there are expected 
to be many more p- than s-wave resonances at the small widths where it was not possible 
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(using means independent of RMT) to differentiate the two parities, and hence an extremely 
large number (see below) of possible "s-wave" sets can be constructed from the observed 
resonances. 

Many different selections were applied to obtain the NDE. For example, in Ref. [2j it is 
stated that "The criterion for inclusion in the NDE is that the individual sequences be in 
general agreement with the GOE." Furthermore, data from many of the included nuclides 
were selected, at least in part, using measures derived from RMT for the GOE. 

atained by the group 



19|), they had "...no 



Data for all but three of the 24 nuclides considered herein were o : 
at Columbia University. According to their publications (e.g., Ref. 
specific tests for s vs p levels, so there may be errors in these assignments." Therefore, 
they relied on theoretical guidance, specifically measures derived from the GOE, to perform 
these separations. For example, for six of the 24 nuclides considered herein, including the 
two having the largest number of resonances, separation of p- from s-wave resonances was 
accomplished 39] by i) assuming all resonances having neutron widths larger than a certain 
(unspecified) size were s wave, ii) calculating the number of s-wave resonances below this size 
by assuming the PTD is correct, and iii) deciding which of the resonances below the threshold 
defined in the first step and needed to achieve the total number calculated in the second 
step, were s-wave by requiring good agreement with four spacing statistics (the Wigner 
nearest-neighboor spacing distribution, p(Sj, Sj+i), A 3 , and the Dyson F test) derived from 
the GOE. Separation of s- from j>-wave resonances for several of the other NDE nuclides 
followed the first two steps described above followed by a Bayesian analysis to decide which 
of the resonances below the threshold defined in the first step to assign to the s-wave set. 
The Bayesian analysis again assumes that the PTD is correct (for both s- and p-w&ve 
resonances) and furthermore requires the average widths for s- and p-wave resonances. The 
latter quantity usually is known only approximately, if at all. Such Bayesian analyses are 
known to be unreliable. For example, several neutron resonances in 64 Zn 15] are known to 
be definitely p wave by their symmetrical shape in transmission (total cross section) data, 
but nevertheless have a Bayesian probability of >99 % of being s wave. 

Given these selection procedures then, it is perhaps understandable that the NDE appears 
to agree well with the spacing statistics examined in Refs. jl, 2| despite that fact that the 
data are neither complete nor pure. Consider, for example, the case of 232 Th, the NDE 
nuclide with the largest number of resonances. The authors of the primary reference from 



14 



which these data were taken 



21 



state that 80% of the s-wave set (i.e. the set included 
in the NDE) were chosen to be s wave because they were too large to be p wave. Hence, 
36 of the 178 232 Th resonances fall below this size and had to be selected as s wave using 
measures derived from the GOE, as noted above. In this same reference, 62 resonances are 
assigned as p wave (and hence by definition are below the s-wave threshold). Therefore, 
there are 98 resonances from which to choose the 36 needed to complete the s-wave set. 
Therefore, the number of possible s-wave sets is astronomically large (~10 27 ), and hence it 
should not be surprising that at least one set can be found to agree with the various RMT 
measures used in Ref. [39j. Similar numbers apply to the other data sets from Columbia. 

Finally, given that a significant fraction of the NDE data were selected using measures 
derived from RMT for the GOE, it seems highly questionable to use the full NDE for any 
test of this same theory. In contrast, the "p-free" results presented herein circumvent this 
problem by eliminating these problematical resonances from the test. 



VII. CONCLUSIONS 



a, 



I have shown that neutron widths in the NDE |2j, when analyzed carefully, reject the PTD 
with hi gh c onfidence. Given that reported deviations from the PTD for individual nuclides 



21 



23 



have occurred on both sides of v = 1, the statistical tests used herein likely 



underestimate the confidence with which the PTD can be rejected for the combined set of 
nuclides in the NDE. Furthermore, I have shown that the NDE is not pure and very likely 
incomplete. Also, measures derived from RMT for the GOE often were used in deciding 
which resonances to include in the NDE. These facts cast very serious doubt on repeated 
claims that the NDE constitutes a striking confirmation of RMT. 

Measurement techniques have improved considerably since the data in the NDE were 
acquired. In particular, new methods for determining resonance spins 4J, |45] hold the 
promise of surmounting the difficulty of separating small s- from large p-w&ve resonances, 
which has been one of the most troublesome barriers to obtaining better data. Unfortunately 
most experimentalists continue to use RMT to correct their data for instrumental effects 
rather than use their data to test RMT. Exceptions to this practice are recent tests 23l. 144 1 
which have revealed significant disagreements between new neutron data and the PTD. 



These new data, together with previously reported disagreements 21 



40h43| . which have 



15 



largely been ignored, are potentially very interesting. Therefore, I urge experimentalists 
obtaining new and improved data to use them, when possible, to test RMT. 
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